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I. INTRODUCTION 

Inspired by the achievement in the study of hadron spectroscopy within the framework of quark potential model [1- 
3] , there have been continuous efforts to apply the quark potential model and the resonating group method (RGM) [4] 
to study nucleon-nucleon, nucleon-meson and meson-meson interactions and scattering [5-20]. Among these efforts, 
the study of kaon-nuclcon (KN) interaction arose a particular interest in the past. Due to the high penetrating 
power of K + meson, it is expected that the study of K + N interaction would provide more information for nuclear 
structures and properties. Since in the K + N interaction, the one-pion exchange is forbidden, the two-pion exchange is 
suppressed [10,11] and there is no annihilation of valence quarks to appear, it is expected that the t-channel one-gluon 
exchange potential (OGEP) plus the harmonic oscillator confinement potential would give a reasonable description 
of the K + N interaction. With this idea, the authors in Ref.[10] calculated the S-wave phase shifts of K + N elastic 
scattering and found that the theoretical results are in quite good agreement with the experimental data in the low- 
energy domain. It is noted that in the calculation, the authors only took the coulomb, spin-spin and Darwin terms 
in the t-channel OGEP without considering the spin-orbital coupling and tensor force terms which contribute to the 
higher partial wave scattering. Subsequently, to investigate P-wave phase shifts, the authors in Ref.[ll] introduced 
into their model the spin-orbital coupling terms originating from the t-channel one gluon exchange and a scalar 
exchange term describing the confinement interaction. As one knows, there are two kinds of spin-orbit coupling terms 
in the t-channel OGEP: the spin-symmetric term and the spin-antisymmetric one. Correspondingly there also exist 
such two terms in the confining potential generated from the scalar exchange. The spin-symmetric terms in the 
two potentials have opposite signs. Therefore, the effect of the spin-orbital coupling is suppressed just as required 
in the study of meson and baryon spectra [2,21]. However, the spin-antisymmetric terms are of the same sign and 
hence would produce large splitting which is in contradiction with the experiment [11]. Such terms are therefore 
dropped out in Ref.[ll]. The calculation in Ref.[ll] showed that if the spin-orbit coupling terms are considered only, 
except for the I = 0, J = h channel phase shift, the sign and magnitude of the other channel phase shifts are well 
reproduced; when the other terms in the t-channel OGEP are taken into account together, there appears a serious 
problem that the 1 = 1 channel phase shifts all become negative, conflicting with the experimental data. Later, the 
S-wave phase shifts of KN scattering are restudied in Ref. [15] by employing the Born order diagrammatic technique. 
In the study, although only the spin-spin coupling term in the t-channel OGEP is considered, the calculated result 
looks fine. Subsequently, the Born approximation was applied to investigate the KN scattering more extensively in 
Ref. [17]. In the investigation, apart from the hyperfine term in the OGEP and the linear scalar confinement, the 
spin-orbital coupling and spin-independent terms in the OGEP are taken into account in the evaluation of P-wave 
and D-wave phase shifts. Nevertheless, the magnitudes of most of the calculated phase shifts are smaller than the 
experimental ones. Particularly, the sign of the theoretical P13 wave phase shift is opposite to the experimental data. 
Recently, the KN phase shifts are recalculated in Ref. [18] within the constituent quark model by numerically solving 
the Hill- Wheeler equation, trying to give a consistent description for the KN interaction and the relevant meson 
and baryon spectra. In the calculation, besides the linear confining potential, the authors used only the Coulomb 
term and the spin-spin interaction term multiplied by a phenomenological coefficient function of Gaussian type. The 
calculated 1 = channel S-wave phase shift is pretty good in comparison with the experiment; whereas for the 1=1 
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channel S-wave phase shift, there appears a big discrepancy between the theoretical result and the experimental one. 
The authors also calculated higher angular momentum phase shifts without including the spin-orbit coupling and 
tensor force terms in the OGEP. Even though the results were considered to be quite reasonable, the calculation is 
not complete theoretically because the spin-orbit coupling and tensor force terms in the OGEP were not taken into 
account. 

From the previous works mentioned above, it is clearly seen that a precise understanding of the KN interaction at 
quark-gluon level still calls for a sophisticated quark potential model which can give a consistently good description 
for not only the KN interaction, but also the KN interaction which has never been investigated in the past. For this 
purpose, it is necessary to incorporate new physical ingredients into the model as suggested in Refs.[ll,18]. In this 
paper, we attempt to investigate the KN and KN interactions and their low-energy elastic phase shifts in a consistent 
way within the quark potential model . The new features of this investigation include: (1) The potential model is 
composed of the t-channel OGEP [1] and the s-channel OGEP [22] as well as a phenomenological confining potential. 
The t-channel OGEP is responsible for the KN interaction, while for the KN interaction where the annihilation and 
creation of a quark-antiquark (qq) pair appear, the s-channel OGEP is necessary to be considered as demonstrated 
in our previous investigations of nn and KK interactions [19,20]. In these investigations, it was shown that the 
s-channel OGEP plays a dominant role for the tttt I — channel S-wave scattering and is necessary to be considered 
for the formation of KK molecular states. As one knows, the two OGEPs are derived from QCD in the nonrelativistic 
approximation of order p 2 /m 2 and contain spin-independent terms such as the Coulomb, velocity-dependent terms and 
spin-dependent terms such as the spin-spin interaction, spin-orbital coupling and tensor force terms. All these terms 
are taken into account in our calculation as should be done in a theoretically consistent treatment. (2) Inclusion of 
the QCD renormalization effect. It is well-known that the OGEPs are derived from the tree diagram approximation of 
the S-matrix elements or the irreducible interaction kernels in the Bethe-Salpeter (B-S) equation. Obviously, to refine 
the potential model, the QCD renormalization effect is necessary to be incorporated into the model. This can be done 
by replacing the QCD coupling constant and quark masses in the OGEP with their effective ones which were derived 
in our previous work in the one loop approximation and a mass-dependent momentum space subtraction [23]. Our 
calculation indicates that the inclusion of the QCD renormalization effect gives an appreciable improvement on the 
theoretical phase shifts, particularly, for the P-wave phase shifts. (3) The contribution from the color octet of the three 
quark cluster (qqq) and the quark-antiquark cluster (qq) to the KN and KN scattering is considered because when 
the kaon (antikaon) and nucleon interact, the color singlet states of the clusters (qqq) and (qq) are possibly polarized. 
This consideration has been justified in the recent studies of meson production and decay phenomena [24,25]. In 
these studies, the color octet of qq cluster plays an essential role in explaining the experimental data. According to 
our calculation, the consideration of color octet can also improve the theoretical results. (4) Nonlocal KN and KN 
effective interaction potentials are derived from the underlying interquark potentials by employing the RGM and used 
to evaluate the phase shifts. We do not use the localized version of the potentials because an inappropriate localization 
would damage the hermiticity of the potentials and induce other unexpected errors. (5) The effect of spin-orbital 
coupling in the t-channel OGEP is necessarily suppressed in the present investigation in order to reproduce the P-wave 
phase shifts. This is consistent with the requirement in the study of hadron spectroscopy. The suppression may be 
achieved by a proper change of the coefficient function of the spin-orbit term in the effective intercluster potentials 
which are derived from the corresponding term in the OGEP. With the considerations mentioned above, we obtain in 
this paper a series of theoretical KN S-wave, P-wave and D-wave phase shifts which are in fairly good agreement with 
the experimental data. In addition, a series of theoretical phase shifts for the KN elastic scattering are predicted. 

The rest of this paper is arranged as follows: Section 2 is used to describe the quark potential model and show how 
to derive the KN and KN effective interaction potentials. Section 3 serves to describe the calculation of the KN and 
KN scattering phase shifts. In the last section, the calculated results are presented and discussions are made. There 
are four appendices. In Appendix A, we show the construction of the color-flavor-spin wave function for the systems 
under consideration. In Appendix B, the effective KN and KN interaction potentials derived in position space are 
listed. In Appendix C, we briefly describe the derivation of the phase shift formula used in our calculation. Appendix 
D is used to make some explanations on the QCD renormalization. 

II. QUARK POTENTIAL MODEL AND EFFECTIVE KN AND ~KN POTENTIALS 

According to the quark model, the KN ( KN ) system may be treated as two quark clusters: the if-cluster (qs) 
( the AT-cluster (qs)) and the Af-cluster (qqq) where q — u or d. The effective KN interaction potential may be 
extracted from the following Schrodinger equation for the interacting q 4 s system by the RGM 

(T + V)^ = E^ (1) 
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where E, T, V and W stand for the total energy, the kinetic energy, the interaction potential and the wave function of 
the system respectively. In the center of mass frame, 
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where T c represents the center of mass kinetic energy, 

^=£(^+^+^) (3) 

i<j=l 

here V^, Vg and V£ denote the t-channel OGEP, the s-channel OGEP and the confining potentials respectively. They 
are separately written in the following. The t-channel OGEP represented in the momentum space is [26] 

t _ iirasCtj r p 2 (m?+m|) 7^\2 , (m.-m,) p , 7\ 

V ij ~ 7= 1 TnT Sm 4 ™ 2 \ q - K) + 2m%m m P ■{1 + k) 

^ {Q-k) 2 'J 'J _^ ^ 

+ (9±Al! + i_ fc).(^,_iL)]_ <l=±£a, ( 4 ) 
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where rriij = mi + nij, a s is the QCD fine structure constant, <Ji are the spin Pauli matrices for i-th particle, C\j is 
the t-channel color matrix defined as 

r t _ r^-^-(^-V)' for ««(««) ^ 

--^--§-, for 

with A a being the Gell-Mann matrix, P, fc and q are the total momentum, the initial state relative momentum, and 
the final state relative momentum of the two interacting particles. 
The s-channel OGEP is [20, 22] 

ys _ ^lEL0i UP I 5(m 2 +m' 2 )-4m' _ 2jf _ 

V V ~ 2mm' ^ 2 ° + °^ ^ Sm 2 m' 2 ^ rrfi m' 2 

-( ( « 2 t TO l 2) + K + -£ ) ^ • ^ +t^(p x g) • - £•) 

~4^(g x j) ' (j - j) - jt^ P ■ ■ ^ (6) 
+ 4^t(P ■ ■ o 3 - k ■ ViP ■ aj +A k ■ (Tik ■ o-j) 
+ J^(P ■ 0*9 ■ o, - q ■ a,P ■ a 3 +4g • o^q ■ aj)] 

where m and m' denote the quark (antiquark) masses before and after annihilation respectively, C-j and F?- are the 
s-channel color and flavor matrices, defined by 

q,- = ^W-Af) 2 (7) 

and 

F^l'iln-Tj +Vr V + + V+Vf + UrUf + U+UJ + fa) (8) 

here Tj are the isospin Pauli matrices for i-th particle, Y t the hypercharge operators, and ( XJf and [/") 
represent the rising and lowering operators of the F-spin ([/-spin) respectively. 

The confining potential, as was done in Ref. [10,11], is taken to be the harmonic oscillator one. In the momentum 
space it is represented as 

^=C| 3 .(2 7 r) 3 M^ 2 V 2 ^ 3 (9 - k) (9) 
where /Xy is the reduced mass of the interaction particles and u> force-strength parameter. 
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Now let us construct the wave function of the KN system from the wave functions of clusters (qs) and (qqq). Since 
there are identical particles between the two clusters, the basis function of the system may be represented as 



$TMsm(pi,P2,P3,P4,P5; P) ^ qqx 

- ^(l--Pu-i , 24-ik)*TMi m (1.2,3,4,5)Ji(pi,W,P3,P4,P6;P) 

where we number the three quarks in the iV-cluster as 1,2,3 and the quark and antiquark in the A'-cluster (or the anti- 
quark and quark in the if-cluster) as 4 and 5, Pji (j = 1, 2, 3) symbolize the interchange operators, \I/ TM i m (l, 2, 3, 4, 5) 

and R(pi,p2,P3,P4,P5', P) represent the color-isospin-spin wave function and the position space wave function respec- 
tively which are constructed from the color-isospin-spin wave functions and the coordinate space wave functions of 
nucleon and kaon. For the KN system, noticing that there is no identical particles between the two clusters (qs) and 
(qqq) , the basis wave function of the system may simply be written as 



$TMsm(pi,P2,P3,P4,P5] P ) = *TMi TO (l, 2, 3, 4, 5)R(pi,p 2 ,P3,P4,P5; P ) (H) 

where f TM i m (l, 2, 3, 4, 5) and R(pi,p2,P3,P4,P5', P) are the color-isospin-spin wave function and the position space 
wave function constructed from the corresponding wave functions of nucleon and antikaon. 

Since the KN system is treated as two clusters, when they interact, each cluster may be in color singlet 1 or in 
color octet 8 as indicated in Refs.[12,20]. Thus, the color-spin-isospin wave function \I> TM i m (l, 2, 3, 4, 5) of the whole 

system may be given by the color singlet part *^^ im (l, 2, 3, 4, 5) or the color octet part i m (^ 2, 3, 4, 5) formed 
by the color singlets or color octets of the two clusters. In principle, we may test a general color structure of system 
under consideration which is given by the following linear combination 

*™i m (l, 2, 3, 4, 5) - a^ M 2, 3, 4, 5) + /Jtfg^Jl, 2, 3, 4, 5) (12) 
where the coefficients a and [5 are required to satisfy 

\a\ 2 + \[3\ 2 = l (13) 

The wave functions ^^ M i (1,2,3, 4, 5) and ^^ M i (1,2,3,4, 5) arc listed in Appendix A. They are constructed by 

antisymmetry of the wave functions of identical particles in nucleon. 

Because we limit our discussion to the interaction in the low-energy regime, it is appropriate to write the position 
space basis function of the KN or KN system in the form 

R(pi,P2,P3,P4,P5;~P) = ^oV(Pl^)^HP2, ~P)^ot ) iP3, ~P)<I>V <J>A, P)^^, P) (14) 

where 4>it\pi, P) and 4>os\pj, P) are the lowest- lying harmonic oscillator states of the A^-cluster and iT-cluster given 
in the momentum space, 

cPit\Pi,~P) = (2V^) 3/2 cxp(-ifo 4 2 pf T *A± p^ ■ p) (15) 

in which P is the vector representing the separation between the centers of mass of the two clusters and parameters 
A± are defined by 

A_=/3 1 = -^— ,X + = (3 2 = ^±^ (16) 
4mi + m,2 4mi + TO2 

here mi denotes the mass of d and u quarks, mi the mass of strange quark. The wave function in Eq.(14) can be 
represented through the cluster coordinates in the form 

R(pi,P2,P3,P4,P5; ~P) = X K (q)X N (k 1 k 2 )T(Q, ~P)Z CM (P) (17) 

where Xk(1 ) and lAiffc^fe) are the internal motion wave functions of the ^(-fQ-cluster (qs) ((qs)) and the iV-cluster 
(qqq) with q and fci,fc2 being the relative momenta in the clusters (qs) ((qs)) and (qqq) respectively, T(Q, p) is the 
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wave function describing the relative motion between the two clusters with Q being the relative momentum of the 

two clusters and Zqm{P ) the wave function for the center-of-mass motion of the whole system in which p is the 
total momentum of the system. According to the RGM, the wave function of the two clusters may be represented in 
the form 

*TMm = J d 3 p$ TM i ms (pi,p 2 ,P3,P4,P5; P)f( (18) 

where $ TM i mj (pi,P2,P3,P4,P5', P) is the basis function defined in Eqs.(lO) and (11) and f(p) is the unknown function 
describing the relative motion of the two clusters. On substituting the above wave function in Eq.(l), according to 

the well-known procedure, one may derive a resonating group equation satisfied by the function f(p). Then, by the 
following transformation 

/(?) = J d 3 RT{p, R)1I(r) (19) 

where 

(20) 



T(~P ~R) - 1 ( 3 fo,3/4 f ,3. <rV 2 fe +<M 



in which b is the harmonic oscillator size parameter, the resonating group equation will be transformed to the following 
Schrodinger equation satisfied by the relative motion of the two clusters 



-^V^*CR) + J d 3 RV(R,R'MR') = eV{R) 



(21) 



where e,/x and ^f(R) are respectively the energy of relative motion, the reduced mass and the Schrodinger-type wave 
function for the two clusters and 

V(R, R') = V\R, R') + V S (R, R') + V C (R, R') (22) 

is the nonlocal KN ( KN ) effective interaction potential in which V t (R, R'), V s (R, R') and V C (R, R') are generated 
by the t-channel OGEP, the s-channcl OGEP and the confining potential. The potential V(R, R') in the Schrodinger 
equation is connected with the potential V{ P , p') appearing in the resonating group equation by the following formula. 

V(R, R') = J d 3 pd 3 P T(R, ~P)V(P, pW, R') (23) 

where V(p,p') is described in Appendix B. To compute the elastic scattering phase shifts, we need to calculate the 
transition matrix between initial and final plane wave functions as follows 

T fi (k,k')= j d 3 Rd 3 R'e- ikR V(R,R')e ik ' R ' (24) 

where k' and k are the KN (KN) relative momenta for the initial and final states respectively. Upon substituting 
Eq.(20) into Eq.(23), it is easy to find 

Tfi (t, k') = I(!§)3/2 e ^U V 2 ) f tfptfp'e-ik-tvtf, ~p )e *'-7' (25 ) 

Z 7TO^ J 

This expression shows that to calculate the transition matrix, we may directly use the potential V(P,p') instead of 
the potential V(R,R'). 
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III. CALCULATION OF PHASE SHIFTS 



The phase shifts of the KN (KN) low-energy elastic scattering are calculated in the Born approximation. As 
argued in Appendix C and demonstrated in the previous literature [14, 15], the Born approximation can reasonably 
describe hadron elastic low energy scattering processes. In this approximation and in the center-of-mass frame, the 
l-th partial wave phase shift is expressed by the following formula whose derivation will be sketched in Appendix C 

Sf J = -2MkT)(\k) (26) 

where 

in which E is the total energy of the KN (KN) system, vhk and mjv are the masses of kaon and nucleon, k =\k\ = \k'\ 

is the magnitude of relative momenta k and k' in the case of elastic scattering and Tj/ l (k) with isospin 7, total 
angular momentum J and orbital angular momentum I is the transition amplitude. This amplitude can generally be 
expressed as 

T I f f\k)= J2 Ci ro8 ^ m J , dn(fe)^(fcO^K* , )^(fc)T/ i (fc,P;m a ,m' a ) (28) 

m ,m / , fj, , fj, f 

A 

where C^ lrn are the Clebsch-Gordan coefficients, Y\ m (k) are the spherical harmonic functions and 

T} t {k,k , ;m s ,m , s ) = (C;/,M;^m 8 | V fi (k,k') \ I, M^^m'^C) (29) 

are the matrix elements of the operator Vfi(k,k') defined in Eq.(25) between the color-spin-isospin wave functions 
| J, Mi; \,m' s \ C) and | /, Mj; |; m s ; C) in which J, Mj and |, m s are the isospin and spin quantum numbers of the 
KN (KN) system respectively and C denotes the color singlet of the whole system. These matrix elements can be 

easily calculated. The explicit expressions of the quantities Tj t (k, k';m s ,m' s ) and Tj\ l (k), we think, are unnecessary 
to be listed in this paper. We only show here numerical results of the theoretical phase shifts in Figs. 1-6 using the 
conventional partial wave notation L/2J- It is noted here that the formula in Eq.(28) is general for evaluating the 
transition amplitude, particularly, in the case that the spin-orbit coupling and tensor force terms are present in the 
nonlocal effective potentials. 

We would like here to discuss the problem of suppression of the effect of spin-orbital coupling. It is a common 
recognition in the study of hadron spectroscopy that the effect of the spin-orbital coupling term in the t-channel OGEP 
ought to be suppressed by the corresponding term in the confining potential [2]. In Ref.[21], one of the authors in this 
paper and his coworkers proposed a qq confining potential which was obtained from a general Lorentz structure of the 
confinement. In the confining potential, there are various terms among which the spin-orbit coupling term is of a sign 
opposite to the corresponding one in the OGEP. In this paper, to avoid the complexity of such a confining potential, 
we alternatively take an effective treatment to achieve the spin-orbital suppression. Looking at the expression of the 

potential V t (p,p r ) shown in Appendix B, one can sec that there is a kind of factorial functions to appear in some 

terms of the potential V*(P, p') which are of the form 

g(x,~p) = [ fL e -^ 2 +xr-^ (30) 
J Airr 

In particular, this function is related to the spin-orbit coupling term in the potential V t (p,p'). The function g(x, p) 
may appropriately be replaced by an interpolating function such that 

^ e (l-7)=£ 

9(x, P)^ 2x (31) 
To obtain the above expression, we have used the approximate expression of the following integral 
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f(x) = f e^dt^xe-^ (32) 
Jo 

As shown in Fig. 5, when we take the parameter 7 = 0.3, the function on the right hand side of Eq.(32) very approaches 
the real value of the integral when x is not too large. However, as shown in Fig. 6, the above value of 7 leads to worse 
P-wave phase shifts. In order to get better P-wave phase shifts, we have to take some larger value of 7 which just 
plays the role of suppressing the effect of spin-orbital coupling. 

Finally, let us discuss the inclusion of QCD renomalization effect. As mentioned in Introduction, the OGEP is 
derived from the S-matrix or the B-S irreducible interaction kernel in the tree diagram approximation. Obviously, 
to improve our calculation, it is natural to consider the correction arising from QCD renormalization. This can be 
done by replacing the QCD coupling constant and quark masses in the OGEP with the effective ones which are 
obtained by solving the renormalization group equations satisfied by the renormalized coupling constant and quark 
masses. This procedure, as proved in Ref.[27] and demonstrated in Appendix D, is equivalent to replacing the free 
wave functions, the free propagators and the bare vertices in the tree diagrams with the exact ones. In the calculation 
of this paper, we employ the effective coupling constant and quark masses given in Ref.[23] which were derived from 
QCD in the one-loop approximation and the mass-dependent momentum space subtraction. These effective quantities 
are suitable to any energy, particularly, to the low energy, unlike the results obtained in the minimal subtraction [28] 
which actually arc applicable only in the large momentum limit. The effective fine structure constant used has the 
expression like this [23] 

= ^ ( 33 ) 

where a° R is a coupling constant and G(A) is a function of variable A which has different expressions given by 
the time-like momentum subtraction (the subtraction performed at time-like renormalization point) and the space- 
like momentum subtraction (the subtraction carried out at the space-like renormalization point). For the time- like 
momentum subtraction, 



G(A) = 11 In A - ^N f [2 + V3tt - A + ( J_ + 1)^-^ ln±(A + \J A 2 - 4)] (34) 

where Nf is the quark flavor number which will be taken to be three in this paper. While, for the space-like momentum 
subtraction, 

G(A) = 11 In A - *N f [£ - 2 - (£ - 1)*™ In I(A + V^TI) 

+V51ni(l + V5)] [ j 



in which A is defined as A = y/q 2 //J? with q being a momentum variable and /1 the fixed scale parameter. The 
expression in Eq.(33) with the function G(A) given either in Eq.(34) or in Eq.(35) will immediately goes over to the 
result given in the minimal subtraction in the large momentum limit. The latter subtraction was performed at the 
space- like renormalization point. It would be noted that in writing the above effective coupling constant, the mass 
difference between different quarks is ignored for simplicity. The behaviors of the effective coupling constants given 
in Eqs.(33)-(35) are described in Fig. 7. From the figures, we see that the effective coupling constants given by the 
time-like and space-like momentum subtractions have different behaviors in the low-energy regime. It is interesting 
to note that the effective coupling constant given in the space-like momentum subtraction is almost the same as given 
in the minimal subtraction in the regime. This situation happens only in the case of ignoring the mass difference 
between different quarks and taking Nf = 3. In other cases, the difference between the results given by the both 
subtractions will be manifest. 

The effective quark mass is represented as 

m R {\) = m R e- s ^ (36) 

where m R is the constant quark mass given at A = 1 which will appropriately be chosen to be the constituent quark 
mass in the quark potential model, 5(A) is a function which also has different expressions for the different subtractions. 
For the time-like momentum subtraction, 

S(X) - ^^{2 + ( Jf - ^) In I 1 - A 2 [} (37) 
While, for the space-like momentum subtraction, 
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S{\) = S 1 (X) + iS 2 (\) 



(38) 



where 



Si(A) = ^[(|- + l)ln(l + A 2 )-4V2] 



(39) 



and 



5 2 (A) = 




[— ln(l + A 2 )-- + l-ln V / 2] 



(40) 



The behaviors of the effective quark masse given by the time-like momentum subtraction and the real part of the 
effective quark mass given in the space- like momentum subtraction are depicted in Fig. 8. The figures show that at low 
energy the effective masses given in the both subtractions are not different so much. For the interaction taking place 
in the t-channel, as explained in Appendix D, the transfer momentum is space-like, while for the interaction in the 
s-channel, the transfer momentum is time-like. Therefore, for the t-channel OGEP, we will use the effective coupling 
constant and quark masses given by the space-like momentum subtraction and in this case, we only adopt the real 
part of the effective masses in our calculation; while, for the s-channel OGEP, the effective coupling constant and the 
quark masses given in the time-like momentum subtraction will be employed. The variable A is usually defined as a 
ratio of the momentum related to the process of quark-gluon interactions. In this paper, as an effective treatment, we 
directly define it as A = k/fj, where k is taken to be the magnitude of the relative momentum of the two interacting 
particles K and N ( or K and N ). 



This section is used to present calculated results for the KN ( KN ) elastic scattering phase shifts, discuss 
adjustments of the theoretical parameters and analyze the effect of color octet and QCD renormalization as well 
as the suppression of the spin-orbital coupling. First, we focus our attention on the KN scattering. The theoretical 
phase shifts of the KN scattering are depicted in Figs. 1-3. In the figures, the solid lines represent the final results 
obtained by considering the contributions arising from the color octet, the QCD renormalization and the spin-orbit 
suppression. To exhibit the effects of the color octet and the QCD renormalization, in the figures, we also show the 
results without considering these effects. Such results are calculated with the same parameters as for the solid lines 
and represented by the dotted and dashed lines respectively in Figs. 1-3. The figures show us that the agreement 
between the final calculated results and the experimental data is good for the phase shifts of all S-waves, Pi 3 wave 
and D13 wave in the low-energy domain, particularly, in the region of the laboratory momentum less than 600 MeV 
to which the nonrelativistic quark potential model is considered to be applicable. For the other P-wave and D-wave 
phase shifts, the agreement is qualitatively reasonable. In obtaining these results, we used the parameters as follows: 
the QCD coupling constant a s — 0.23, the constituent quark masses m u = = 350 MeV and m s = 550 MeV, 
the size parameter of harmonic oscillator b = 0.255 fm, the force strength of confinement u> = 0.2 GeV , the color 
combination coefficient a = 0.915, the scale parameter of QCD renormalization /i = 0.195 GeV and the parameter 
of spin-orbital suppression 7 = 0.45. These parameters are adjusted to give a better fit to the KN elastic scattering 
experimental data, mainly to the S-wave phase shifts because the KN elastic scattering data are available and rather 
sufficient [29,30]. In comparison with the previous results given in Refs.[I0, 11, 17, 18], our calculation achieves a 
considerable improvement on the theoretical phase shifts for all the partial waves not only in the magnitude, but 
also in the sign. Especially, for the P13 wave phase shift, it now gets a right sign in our calculation, opposite to the 
previous result which was given a wrong sign [17]. 

Now let us analyze the effects of the color octet, the QCD renormalization and the suppression of the spin-orbital 
coupling. In adjusting the theoretical parameters, we found that the calculated results are sensitive to the parameters 
a and b. Any small change would cause noticeable influence on calculated values. For instance, when we let a increase, 
the absolute values of S-wave phase shifts decrease rather fast. In particular, when the a tends to unity, i.e. the color 
octet is absent, as denoted by the dotted lines in Fig.l, we obtain the S-wave phase shifts similar to those given in 
Rcf. [17]. In this case, certainly, we may give a better fit of the calculated result to the experimental one by adjusting 
the parameter b and others, but, we failed to simultaneously get a good result for another S-wave phase shift, as was 
demonstrated previously in Ref.[18]. Only when the color octet is considered, it is possible to get consistently good 
results for the both S-wave phase shifts as denoted by the solid lines in Fig.l. This suggests that the introduction of 
the color octet is necessary in our calculation. From Figs. 2 and 3, we also see that the inclusion of the color octet 
gives a appreciable effect on the P and D wave phase shifts. 



IV. RESULTS AND DISCUSSIONS 
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In this paper, the two kinds of spin-orbital terms in the t-channel OGEP: the spin-symmetric term and the spin- 
antisymmetric one are all taken into account. In this case, we still encountered the puzzling problem for the P-wave 
phase shifts as was revealed originally in Ref.[ll] and mentioned in the Introduction. If the QCD renormalization 
effect is not considered, except for the Pn and P13 wave phase shifts, it seems to be able to reproduce the phase shifts 
for the other phase shifts by readjusting the parameters involved. However, as exhibited in Fig. 2, for the Pn and Pi 3 
waves, it is impossible to get a satisfactory result of their phase shifts. Particularly, for the Pi 3 wave, its phase shift is 
always of a wrong sign as the previous result given in Ref.(17). This problem can only be resolved by taking the QCD 
renormalization effect into account in our calculation. In this way, the Pn and P13 wave phase shifts are accessible 
to the experimental values. As shown in Fig. 3, the QCD renormalization effect gives an essential improvement on 
the D-wave phase shifts as well. Particularly, it renders the D i3 and D i5 wave phase shifts to have the right signs. 
But, this effect is not noteworthy for the S-wave phase shifts. This explains why the previous investigations could 
give some rather reasonable results for the S-wave phase shift. However, as shown in our calculation, in order to get 
the desirable phase shifts for all partial waves, it is necessary to incorporate the QCD renormalization effect into the 
model used. In addition, to achieve such results, as mentioned in Sect. 3, the spin-orbital coupling effect is necessary 
to be suppressed. The necessity of the suppression is separately illustrated in Fig. 6 for the P-wave phase shifts only. 
This is because the spin-orbit term in the effective potential gives no contribution to the S-wave scattering and it 
mainly affect the P-wave scattering. Fig. 6 indicates that when the parameter 7 is taken to be the value 7 = 0.3 

which makes the function g(x, P ) reach its real values, the P01 wave and the Pn wave phase shifts are far from the 
experimental ones, but, when the 7 is getting larger, the absolute values of the P01 wave and the Pn wave phase 
shifts become to be comparable with the experimental results. For the other P-waves, the phase shifts evaluated at 
7 = 0.45 are also better than those given at 7 = 0.3. 

Let us turn to the KN scattering. At present, the low energy elastic and inelastic experimental data for the KN 
scattering are insufficient [31-33]. Therefore, the detailed partial wave analysis for the scattering is almost absent. 
But it is a common conclusion suggested in the previous investigations that the K~p interaction is strongly attractive 
[34]. For the KN scattering in the I = lchanncl, i.e., for the K~ — neutron scattering, there is almost no available 
data and different theoretical models give different predictions. In view of this situation, our calculated results for 
the KN phase shifts can only be viewed as a theoretical prediction [35,36]. It would be noted that unlike the KN 
interaction for which the exchanged part of the effective potential generated from the t-channel OGEP is dominant, 
for the KN interaction, there is no such an exchanged potential. Instead, the direct part of the s-channel OGEP 
plays an essential role and leads to an attractive interaction as seen from the positiveness of the phase shifts plotted 
in Fig. 4. In the figure, the prediction for the P and D wave phase shifts are simultaneously given as well, all the 
phase shifts in the figure are presented in the momentum region less than 200 MeV where any resonance could not 
appear. Here we take the K~p S-wave phase shift given in this paper as an example to estimate the reasonability of 
our calculation. From the relation y/ir | a e ^ |~| f + | [37] where a e ^j is the effective scattering length and f + is the 

i-th wave scattering amplitude and the formulas given in (C.2) and (C.15), we have 5j J ~ kRe{a l e ^). By this relation 
and the calculated S-wave phase shift, it is found that the real part of the scattering length is above 0.58 which just 
lies in the range shown in Rcfs. [33,35]. In addition, we note that our calculation predicts a weaker attraction for the 
K~ — neutron interaction which is different from previous result (see nucl-th/0004021). 

At last, it should be emphasized that the quark potential model used in this paper was established in the nonrel- 
ativistic approximation of order p 2 /m 2 , therefore, the calculated results are only valid for the KN and KN elastic 
scattering in the low energy domain. The model used can not give a complete description for the inelastic scattering 
and the production of resonances which would appear for some higher partial waves in the higher energy regime. 
To explore the KN and KN inelastic scattering in the higher energy regime, it is necessary to apply a relativistic 
approach or a nonperturbative theory. Anyway, the investigation based on the constituent quark model is meaningful 
as it is not only helpful to understand KN and KN interactions from the underlying dynamics, but also provides a 
firm basis of studying five-quark bound states. 
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VI. APPENDIX A: THE COLOR-FLAVOR-SPIN WAVE FUNCTIONS 



In general, the color singlet color state of the five quark cluster (q 4 s) or (q 3 qs) may be built up by the color 
singlets of the TV-cluster (qqq) and ^-cluster (qs) (or the TV-cluster (qs)) or the color octets of the two subclusters. 
Correspondingly, for the five quark cluster, there are two classes of color-flavor-spin wave functions denoted by 

^™Li (1,2,3,4,5) and *™Li (1,2,3,4,5) which are color singlets as a whole, but associated respectively with the 

_/ Ivl ^ 777 J_ Ivl ^ 777 

color singlets and the color octets of the two subclusters. In the function ^J^. im (l, 2, 3, 4, 5), the color-flavor-spin 
(CFS) wave function $^1^ (1,2, 3) at for the TV-cluster which is totally antisymmetric ( of the symmetry denoted 

by the Young diagram [l 3 ] c /s ) is constructed from the C-G coupling of [l 3 ]c x [3]fs where [l 3 ]c and [3]_fs are the 
Young diagrams denoting the antisymmetric color singlet and the symmetric flavor-spin states respectively. In the 
function ^^j f im (l, 2, 3, 4, 5), the antisymmetric CFS wave function ^±) M \ m (1,2, 3)jv for the TV-cluster is given by 

the C-G coupling of [21] c x [21]fs where [21] c and [21]fs represent the color octet state and the flavor-spin state of 
mixed symmetry respectively. The explicit expressions of the wave functions mentioned above can easily be written 
out by the familiar method given in the group theory, as displayed in the following. 
The first class of the CFS wave function in Eq.(12) for the whole system is 

M 1 M 2 

where ^ij^ i m (l> 2, 3)jv , as mentioned before, is the CFS wave function for the TV-cluster and *i^ 200 (4, 5)k is the 
CFS wave function for the if -cluster. They are represented separately as 

^ 1 i ms (l I 2,3) Ar =e c °(l,2,3) X «, 



*iiL i™. (L 2 > 3 V = £c°(l, 2, 3)xil , (1, 2, 3) (A.2) 



where 



£(1, 2, 3) = ^=e abc q a (l)q b (2)q c (3) (A.3) 

represents the color singlet wave function of the TV-cluster and 

xf Ml , ms (1, 2, 3) = -L [ X \ Mi (1, 2, 3)^ ms (1, 2, 3) + X \ Ml (1, 2, 3)^ (1, 2, 3)] (A.4) 

is the isospin-spin wave function of the TV-cluster in which the isospin wave functions Xi Ml (1, 2, 3) and Xi Ml (1, 2, 3) 
and the spin wave functions f1 m (1, 2, 3) and <y9^ m (1, 2, 3) are expressed as follows 

X| Ml (1,2,3) = E ^ m3 Cj™ i , m2X , mi (l) X , m2 (2) Xim3 (3) 

2 m,m3mi,m2 2 * z 



X\ M1 (1,2,3) = E 4TM 1 ^° mi > 2 X> 1 (l)X| ro2 (2)X| m3 (3) 

z mi ,m2 ,m.3 2 j. j. 

^ ros (l,2,3) = E ^^^,^(1)^,^(2)^^3(3) 

2 m, 7713, mi ,771 2 

4 ms (l,2,3) = E ^^^^(^^^(^^(S) 



(A.5) 



The CFS wave function of the TV- cluster is 



*|moo( 4 < = C o( 4 ' 5 )xim(4, 5 W4, 5) (A.6) 

where Co(4, 5), XiM(4,5) and yoo(4, 5) are the color, isospin and spin wave functions, respectively. Since there is no 
identical particles in the cluster, these wave functions are of the forms 

C (4,5) = -^(4)<f(5) (A.7) 

and 
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Xi M (4,5)= E C|^ i00 xi mi (4)xoo(5) 

mi m 2 2 (A 8) 

^oo(4,5)= E GS° mi4ma Vi mi (4)v ima (5) 

mi m2 

For the second class of the CFS wave function in Eq.(12), it can be represented as 

*S,|JM,3,4,5)= £ £ C^ M2 *g iim (l, 2, 3) N *( 2 ); oo ( 4 , 5 ) x (A.9) 

MiM 2 c 

where * ( x 2 ^ x (1,2,3V and *r?L,(4,5V are the second class of CFS wave functions for the iV-cluster and the 
if-cluster respectively. Their expressions are shown in the following. 

^iSi™.^ 2 ' 3 )^ = ^[^a^'^fSi^/l.a.S)-^ (l,2,3) X ^ iTO .(l,2,3)] (A.10) 

where #(1,2,3) and #(1,2, 3) are the color octet wave functions given respectively by the Young- Tableau [211] and 
the Young- Tableau [121] and x^ilt i (1)2, 3) and x^ilf i (1j 2, 3) are the corresponding isospin-spin wave functions. 
Their expressions are 



#(1,2,3) = |e« 6 [g°(l)g i (2)^(3) + g°(2)g i (l)^(3)] 

#(1,2,3) - 2j 5 c«6[g°(l)« < (2)^(3) -g°(2)« i (iy(3) - 2 g»(3) ? i (l)^"(2)] 

xS£ , ros (1, 2, 3) = ^ [ X | Mi (1,2, 3)^ m> (1,2,3)- X \ Ml (1, 2, 3)^ (1, 2, 3)] (A-") 
xf^ im , (1, 2, 3) = -^[X| M1 (1, 2, 3 V| ma (1, 2, 3) + X \ Ml (1, 2, 3)^ m< (1,2, 3)] 

The second class of the CFS wave function for the if (or K )-cluster is as follows 

*imoo( 4 ' 5 )- =Ca(4,5)xiAf(4, 5Voo(4,5) (A.12) 

where 

C a b (4, 5) = q b (A)q a (5) - ±6 b a q c (4)q c (5) (A.13) 
is the color octet for the K(K) cluster and the other two functions Xam(4> 5), </?oo(4, 5) are the same as in (A. 8). 



VII. APPENDIX B: THE EFFECTIVE KN AND KN INTERACTION POTENTIALS 

In this appendix, we show the nonlocal effective interaction potentials of the KN and KN systems which are 
derived from the intcrquark potentials and the RGM. 

The KN nonlocal effective potential V t (P,p') which is derived from the t-channel OGEP written in Eq.(4) is divided 

into two parts: the direct part V t D (P,p') and the exchanged part V t ex ( P , p')\ 

V t (t,p')=V t D (p,p')-VL(p,p') (B.l) 

where 

vr(p,p') = vr(p,p') 14 + vr(p,p , ) 24 + vr(p,p') 34 (b.2) 

here the superscript ab = 14, 24 or 34 designate which pair of quarks interchange. Each part of the potential contains 
several terms as shown in the following 



v t D {P,p')^vg{P,p\+v£(P,pi) + vi{P,p>) (B3) 
+v 1 d A^,p') + v 2 d 4 (p,p') + v 3 d 4 (p,p') 



n 



vr ( p , pT = ( p , p'Y h + v 2 7 ( p , pT + Viz ( p , P 'Y b + v& ( p , P 'Y b + v 2 f ( p , P ') ab (B 4) 

+v 1 f(p, P ') ab + v 2 f(p, p') ab + Vfg(p, P 'Y h + v 4 f{p, p'Y b + vg(p, P 'Y b 

the subscript in each term on the right hand sides (RHS) of (B.3) and (B.4) marks the two interacting quarks: one in 

the TV-cluster and another in the ^-cluster. The terms V®{p,p') and V^ x (p , p') ab are derived by the RGM in such 
a way 

v£(p,7') = J tl^^(R(pi,P2,Ps,P4,P 5 -, ~p) I I rAmMip')) (B.5) 



and 



5 — *■ 

v^(p,p'Y b = J f[ {Zf ( r (p^p^p^ ~p) I v ii p ^ I i?(?i,5,5,i,i;?)> , (B.6) 

where the quark potential V^* was denoted in Eq.(4) and the position space wave function was given in Eq.(14). 
First we describe the ten terms on the RHS of (B.4). By introducing the following functions: 





-xr 2 


fir • P 




fL(P,P')ex = e~ 


p V 

2b 2 


-£-^<?-?> a 




fk(P,p')ex = er 


p y 

2b 2 


-^(p+?) 2 -^(p- 


-p') 2 


ff 5 {P,p')ex=e- 


p p' 


-^p 2 -ff(pV) 2 




fl4(P,p')ex=e- 


p p' 

2b' 2 


-^(p-?) 2 -^(p- 


V) 2 


f2 3 {P>p')ex=e- 


P V 
2b 2 


-ffHp-?) 2 





(B.7) 



the exchanged terms of the potential V t ex (P , p') lA can be written as: 

/\14 _ iira s Cl 



V 2 f(P,pT = ^^fL(P,p')ex{9(^ P) ^(1+ o 2 ■ a A ) -^[fc 

P - | a 2 • a 4 + 1){~P x ?) • 2 - a 4 ) (B ' 8) 

(l+ 7 2 ) p ■ a 2 p ■ a 4 ]g(^,p)} 



v&(r,K) 14 = ^0frfLCp,p'U{g{Tw, ~p) - 4^(i+ ^ • * 4 ) - 

P pT - 4 ^3 • ^4 +»&(1 + 7)(P x p>) ■ (a 3 a A ) ( B - 9 ) 

(I + 7 2 ) P • ^sP • <?4}g( 2 j?,P)} 

e.x(~X _/\14 47ra s C 



14 



(P,P') 14 = ^f^A t 4(P,P')e ;C {.9(^, P - P') - 4^(!+ ^1 • ^) 



-i^t-^i^^- P'f + L£ ¥ )1 -9-1^.^4+^(1+ 7 )(p x p') • + a 4 ) (B.10) 
+T^(1 +7 2 )(P - ?)■ *i (P - ?)■ ^M^, P - ?)} 

V?2*(P,pV = S#^4(?, P)e,{5(^,?) - 4^(1+ -1 • -2) - J^i^ 

-M 2 ^ P' -h ~P) 2 J~ 2 ^ ' ^ +i^r(l +7)(? x p) • - a 2 ) ( BA1 ) 

+w(1 + 7 2 )p'-^ip'- H.g(^,p')} 
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V ex Co n'\ 1A — 47ra ° C 45 ft I J 7j\ I„(a 2 n i\ _ / ml+ml a 4 -a 5 \ 1 [6 

v Ah\>iP) ~ ( 7Tb 2/ ai - ) 3/2.ll5\P >^ )exW(-p-,P ) \ 8m 2 m 2 i 4mira 2 > 4( mi +m 2 ) 2 b^ L Ql 

(aip'-/3ip'+ftP) 2 a 2 ft , 

1 \g n . 02(01-02)^1 rt ./ Q2(ai-/?l) 2 -a2/3lfe / 
' 4m im2 ^ l b °^ "j^p P ■ ^2 P ~a 2 CT 4 • CT 5 



(B.14) 



? P)^- i *i • *3 + 7)(? x P) • (a! - a 3 ) ( B - 12 ) 
+ T i E t(l+~ f 2 )p : '-* 1 p'-aM^,P')} 

V^r P T = ^^(P,?).^ - ^ - ^ + ^(1 W)CP - ?)»] (B.13) 

T/ ea 7^7 _/\14 _ 47ra 3 C; 5 ft ^ a r / o 2 7^ _ ^ ™i+™2 , "1^5 1 1 [ 6 

v VS\l J ,P) - ( 7Tb 2/ ai }3/2.ll5(r , P )exW(jT, V ) \ 8m 2 m 2 -I- 4mim 2 > 4(mi+m 2 ) J 6 2 Lqi 

_ (oi^p+^T _ o^ i(1 _ 7)( p x p } . [(1 + ^ } - _ (1 + ^ } a B ]]g(^J) 
-3^[6«2 - p • ? _ £fe^£^ ; 2 _ a2 ^ . ^ 5 

V25 (P,P') 14 = J^W^fT2 fkiP ^ P')ex{9{^r P + P') - (TmfTf + 4mim 2 ) 

- 4( mi+ U)^ ^ - - + jfe^ 1 " 7)(P x ?) ■ [(1 + a 2 

-(i + ^) * B ]]s(f* , p + ?)]Mt> p + ?) - 5^[ 6 «2 + f(«i - fr 

-2a 2 /3 2 )(p - ?) 2 - a 2 o 2 ■ a 5 + 7 )(p^x p') ■ [(2 + ^) ^ 2 

+ (2 + 21) a 5 ] - if (1 - 7 2 )(p + h °2 (P + ?)• ? 5 Mf , P + ?)} 

- ^w ^ - - + ^^fe^ 1 - 7)(P x ?) • [(1 + jj*) a 3 

-(2 + ? 5 ]Mf , P + ?) - i^[6^ + f («i - f - 2«2/?2)(P - pT ( BA6 ) 
-ai ^3 • ^5 + 7)(P x p') • [(2 + ^ 3 +( 2 + a) a 5 ] 

-#(1 - 7 2 )(P + ?)■ ?3 (P + ?)■ * bMS*. P + ?)} 



(B.15) 



'^ lP)2 ~ gfr<(l - 7)(P' x P) • [(1 + a 4 -(1 + ?5]]0, P') (B17) 



in which a\ — TOi ™ 1 TO2 , Q-i = m "+ m2 ■ The other exchanged potential V t ex (p , p') 24 and V t ex (p , p') 34 can directly be 
written out by changing the superscripts 14 to 24 and 34. 

The terms of the direct part of the potential in (B.3) are shown below 

V 2 d a CpJ) = p^/^fPyiDls^,^ ?) + 3^(1- ^ 2 • a 4 ) - 4^ [6- i ^2 • ^4 

-^3^ P 2 +^(! + 7)(P x ?) • + - ^(P - ?) 2 - 7)(P x ?) ( R18 ) 

•(^2 - ^4) + 1^(1 + 7 2 )(P + P')' ^2 (P + P')- ^4]5(2^, P + P')} 

where 

/| 4 (P , = e-^-^-T^V) 2 (B.19) 
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The terms V^(p,p') and V^(p,p') in (B.3) have the same form as shown above except for the subscripts 24 being 
changed to 14 and 34. 

The term V£(P,p') in (B.3) is of the form 

V8(t> ?) = T$7£fa&(t> P')n{9(^J + 7')- (£gf + & 

-(1 + ^) ^LMfJ, P + ?) - 4^^[6«2 + f*(ai - £ - 2 «2/?2)(P - ?) 2 ( B - 2 °) 
-a 2 ^2 • * 8 + 1){~P x ?) • [(2 + ^) a 2 +(2 + 2L) a 5 ] 

-#(1 - 7 2 )(P + ?)• ^2 (P + ?)■ ? 5 ] 5 (f , p + ?)} 



where 



/l 5 (P,P% = e-^l'VMK'V)' (B.21) 



The remaining two terms in (B.3) are of the same form as given above except for the subscripts 25 being replaced by 
15 and 35. 

Let us turn to the KN interaction potential. For the KN interaction, there is only a direct part of the potential 
coming from the t-channel OGEP as represented in (B.3), (B.18) and (B.20) because there are no identical particles 
between the iV-cluster (qqq) and the if -cluster (qs). In addition, the nonlocal effective potential derived from the 
s-channel OGEP plays an essential role in the KN interaction. This potential can be written as 

V s (~P, P') = V^ d (p, p') + V 2 S 4 "V, p') + V^ d (p, p') (B.22) 

here Vf^ d (P, p') denotes the direct term of the potential generated from the interaction of the quark 1 and the 
antiquark 4 

V?4-V,P0 = 4 ^ffi /»(^){(3+ *i • ?4) - ^[3 - ^iH(P - p'f] 

-5^(2+ ?i • - 4 )(3 - igs£(P x p>) ■ (a, a 4 ) (B.23) 



+ ~ J^\P ■ ■ ^ + ft ■ ■ ^4]} 



where 



^ ^ 3g 2 -l/2 - ,_ 3/3 2 +l/2 / N 

fL{P^)n = e 2b ? 4b ? (B.24) 

The other two terms V£± d {P -,p') and V^4- d (p,p') can be written out from the above expression by the substitution of 
the subscripts 24 and 34 for 14. 

The effective KN (KN) potential derived from the interquark harmonic oscillator confining potential and quark 
interchanges arc of simple expressions. They are represented as 



V c ex (P, P') = -12& 2 ^ 2 {Cf M 24 + C S 34 M 34 + Cl 2 p 12 + C 13 M 13 + C s 14 /!14 



P V 3<3 2 



( p -p') 2 



(B.25) 



where the p,ij denotes the reduced mass of interacting quarks i and j and the color factors are the same as the ones 
appearing in the t-channel potentials. 

Apart from the potentials listed above, there are additional terms in the KN and KN potentials occurring in the 
resonating group equation which arise from the kinetic term and the normalization term in the equation due to the 
effect of quark rearrangement. They are written as follows 

= WcVH&fr ' ( :f^ +p,)2 ] + - i 1 (B.26) 

J 3_[ 1 P_l 4. J_[3a2. _ ajp' 1-1 fexfZ n '\ 

^2n 2 LIbZ 1446 4 J ^ 2n 3 \-2 b 2 46 4 \5Jt\ i 'P) 
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where 



frCp.P 1 ) = e-^-^ (p - p ' )2 (B.27) 



/" = — -, ; .Mi = ^-,M2 = -5-.M3 = 1 ,7i = 3/(1/3 + ai). B.28 

4mi + m 2 2 3 mi + m 2 

and 

7V e:E (p,p') = -3^re"^ i^( p -") (B.29) 

here E r is the relative energy of two clusters. 

The color-spin-isospin matrix elements of the above potentials are easily evaluated by using the CFS wave functions 
given in Appendix A. 



VIII. APPENDIX C: DERIVATION OF PHASE SHIFT FORMULA 



In this appendix, we briefly describe the derivation of the formula used to compute the phase shifts with a comment 
on the Born approximation. One of the authors of this paper, in his previous work on the rclativistic Pauli-Schrodingcr 
equation for two-body scattering states, which was proved to be equivalent to the corresponding Bethe Salpeter 
equation [38] , proved that in the relativistic case, the outgoing state wave function of the system under consideration 
may be written, in the position space, as 

yj(r) = ^(r) + f(n k )^ (C.l) 

where ip (f ) is the wave function for free particles and f(flk) the probability amplitude of the outgoing spherical 
wave . For a boson system, the wave function in (C.l) is scalar, while, for a fermion system, the wave function is 
represented in the Pauli spinor space. In the unequal mass case, the amplitude /(fifc) is related to the transition 
amplitude Tfi in such a fashion 

m k ) = -^T fi (C.2) 

where 

M(E) = Ei - {m « E ; m " )2 (C.3) 

with E, mx and being the total energy of the system, the kaon mass and nucleon mass and Tfi is the exact 
transition amplitude. The amplitude is defined by 

T fi = I V I Vi) (C4) 

where ip® is the plane wave function of final state, V stands for the interaction Hamiltonian operator and ipi is the 
exact initial wave function which is determined by the following equation 

1>i = cp* + G°ViPi (C.5) 

in which ip ° is the plane wave function of initial state and 

G° = ] ,e^0+ (C.6) 

E-Ho + ie' ' 

is the Green's function with Hq being the free Hamiltonian. The solution of equation (C.5) can formally be represented 
as 

n=0 
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In the lowest order Born approximation, ipi = Lp\ and, correspondingly, the transition amplitude in (C.4) takes the 
form as given in Eq.(24). 

Now we proceed to derive the formula written in Eq.(26). Let us expand the wave functions in (C.l) in partial 
waves 

^ = ^^P<(cos0) (C.8) 
i 

/(0) = ^/^(cos0) (C.9) 
i 

and 

V- = e lk ^ = ^(21 + l)i l ji(kr)Pi(coB0) (CIO) 
i 

where Pi(cos9) is the Legendre function of rank I. With these expansions, noticing the asymptotic behaviors of the 
spherical Bessel function ji (kr) and the function ipi (r) 

m^y^Hkr- 1 ^) (C.11) 

Mr^yMkr-^ + St) (C.12) 
One may find a result from (C.l) such that 

/* = ^V i5l -l) (C.13) 

On the other hand, the transition amplitude may also be expanded in partial waves 

T/i = 471-^(2/ + l)T,f}(cos0) (C.14) 
/ 

where 

1 f 1 

Ti = —J dxPi(x)T fi (C.15) 

Combining (C.2), (C.9), (C.13) and (C.14), it is easy to derive the following relation 

e 2i5, = 1 _ 4iM(E)kTi (C.16) 
When (C.7) is substituted into (C.4), we can write 

oo 

T fi = J2 T fi (C- 17 ) 

n=0 

where 

T£ = (<p° f | (VG°) n V | 0°) (C.18) 

On inserting (C.17) into (C.15), we have 

oo 

Ti = Y, T ? ( c - 19 ) 

n=0 

where 



1G 



(C.20) 



Upon substituting (C.19) and the Taylor expansion of e 2tSl into (C.16), one may find 

oo oo 

]T ^ {2iS) n = -%AM{E)k J2 (C21) 

n— 1 n— 

From the above equality, it is clearly seen that there exists an one-to-one correspondence between the terms of the 
same order in the both series. When only the first term in each series is considered, we obtain 

Si = -2MkTi(k) (C.22) 

which is proportional to the interaction Hamiltonian V. This just is the formula for the phase shift which is given in 
the so-called Born approximation. If the phase shift is really proportional to the interaction Hamiltonian V as in the 
Born approximation, from the corresponding higher order terms of the both series in (C.21), we can write 

« -^M^T? (C.23) 

which is proportional to V n+1 and means that the n-th terms in the both series in (C.21) are approximately equal to 
each other as if (C.23) is approximately given by taking the equality in (C.22) to the (n+ l)-th power. If the relation 
in (C.23) holds, we see, the formula in (C.22) appears to be a good approximation. 



IX. APPENDIX D: NOTES ON QCD RENORMALIZATION 

To help understanding of the renormalization formulas written in Sect. 3, in this appendix, we give some explanations 
of the QCD renormalization. It is well-known that the t-channel OGEP in Eq.(4) and the s-channel OGEP in Eq.(6) are 
usually derived from the lowest-order S-matrix elements given by the tree Feynman diagrams representing respectively 
the quark-quark(or quark-antiquark) scattering and quark-antiquark annihilation processes in the nonrelativistic 
approximation of the order p 2 /m 2 . As an example, we write the lowest order S-matrix element for two-quark scattering 
as follows 

4? = « (0 Hpi)^r a U (°)(q 1 )z4 J ab (fc)^ )(p 2 ) i5 7 ,y T fc U ( )( g2 ) (D.l) 

where u(°\p) is the free quark wave function, u^°\pi) is its Dirac conjugate (here the spin indices of the spinor 

functions are suppressed for simplicity), iD^} ab {k) is the gluon free propagator with k = pi—qi = q 2 — P2 and ig^T a 
is the bare vertex with g being the coupling constant and T a the color matrix. Correspondingly, the exact S-matrix 
element given by the one-gluon exchange interaction is represented as 

S fl = ^p 1 )T^(p 1 ,q 1 )u(q 1 ) l D^(k)u(p 2 )T^(p 2 ,q 2 )u(q 2 ) (D.2) 

where u(p), iD°^ v {k) and T afl (p, q) denote the full (or say, dressed) quark wave function, gluon propagator and quark- 
gluon vertex respectively in which all higher order perturbative corrections are included. According to the well-known 
renormalization relations 

u(p) = ^fZ^u R {p),u{p) = sfZ^u R {p), , » 

Df v {k) = z.Df^k), r-"(p, q) = z r r%( P , q) V- 6 > 

where the subscript R marks the rcnormalizcd quantities, \/~Z 2 , Z-$ and Zr = Z 2 ~ 1 Z 3 2 are the renormalization 
constants for the quark wave function, the gluon propagator and the quark-gluon vertex respectively, one can get 
from (D.2) that 

S fi = u R {pi)T a £{p u qi)u R { qi )iD^ v {k)u R {p 2 )Y R v {p 2 , q 2 )u R (q 2 ) (D.4) 

As shown in Ref (27), the renormalizcd quantities in the above can be determined by their renormalization group 
equations and their renormalization boundary conditions. The results given by solving the renormalization group 
equations are 
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u R {p)=e^i^ {X) uf{p) 

u R (p)=e?fi^ W u%\p) (D5) 

where (p) and (p) arc of the same forms as the free wave functions except that the quark mass in them becomes 
the effective (running) one, iD R ^ / b (k) formally is the same as the free propagator but the gauge parameter in it is 
replaced by the effective one, <7_r(A) in the bare vertex ig R (\)j^T a is the effective coupling constant, ^72(A), 73(A) 
and 7r(A) are the anomalous dimensions defined by 



= X 7x ln V^W, 73(A) = A A lnZ 3 (A), 
7 r(A) = AAi n z r (A) = -72(A) - ± 73 (A) 



On inserting (D.5) into (D.4), we see that the anomalous dimensions are all cancelled out with each other. As a result, 
we have 

Sfi = 4 0) ( Pl )* 5 «(A)7^r^W( g 0^^^ fc )4 0) fe)^(A)7^^W( g2 ) (D.7) 

This S-matrix element is completely the same as the one shown in (D.l) except that the quark mass, the gauge 
parameter and the coupling constant are replaced by the running ones. In the nonrelativistic approximation of order 
p 2 /m 2 7 one may derive a OGEP from (D.7) which in the Feynman gauge is just as that written in Eq.(4) with 
the coupling constant and quark mass being the effective ones. These effective quantities precisely represent the 
renornalization effect and at one-loop level, they are of the forms as given in Sect. 3. For the s-channcl OGEP, the 
discussion is completely the same. 

Next, we would like to address the renormalization point. In the ordinary QCD renormalization performed in the 
minimal subtraction scheme, which is suitable in the large momentum limit because only in this limit the quark mass 
can be set to be zero, the renormalization point is chosen to be space-like. This choice is suitable to the scattering 
process because in this case the transfer momentum, i.e. the variable of the gluon propagator in (D.2) is space-like. 
This can be seen from the following derivation, 



k 2 = (pi - qi) 2 = 2m 2 - 2sjp? + m 2 ^qt 2 + m 2 + 2 \pt\ \q?\ cos9 (D.8) 
in the high energy limit, we can set m w 0, therefore 

k 2 w -2\pt\ \qt\ (1 - cos(9) < (D.9) 



In the low energy domain, since 



(D.8) can be approximated as 



t 2 



+ m 2 Kim+ J_ (d.io) 



k 2 « -(pt-qt) 2 <0 (D.ll) 



So, for the t-channel OGEP, it is suitable to use the effective coupling constant and quark mass given by the subtraction 
performed at space-like renormalization point . While, for the annihilation process, the momentum in the gluon 
propagator is time- like because in this case k = pi + p 2 = qi + <?2 , 



k 2 = (pi +P2) 2 = 2m 2 + +2\Jpf + m 2 ^ 2 + m 2 - 2 \pt\ |p2|cos<9 (D.12) 
In the large momentum limit, we can write 

k 2 w 2\pt\ \pZ\ (1 -cos9) > (D.13) 

in the low energy regime, we have 

k 2 w {pi - pZ) 2 > (D.14) 

Therefore, for the s-channel OGEP, it is appropriate to use the effective quantities given by the subtraction carried 
out at time-like renormalization point. 
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XI. FIGURE CAPTIONS 

Fig.l: The theoretical KN S-wave phase shifts in the 1 = and 1 channels. The solid lines represent the phase shifts 
with considering the effects of the color octet, the QCD renormalization and the spin-orbit suppression. The dotted 
lines denotes the result without considering the color octet and the dashed line shows the result without considering 
the QCD renormalization. The experimental phase shifts [29, 30] are shown by black squares with error bars. 

Fig. 2: The theoretical KN P-wave phase shifts in the 1 = and 1 channels. The solid lines represent the phase shifts 
with considering the effects of the color octet, the QCD renormalization and the spin-orbit suppression. The dotted 
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and dashed lines denote the results without considering the color octet and the QCD renormalization respectively. 
The experimental phase shifts [29, 30] are shown by black squares with error bars. 

Fig. 3: The theoretical KN D-wave phase shifts in the I — and 1 channels. The solid lines represent the phase shifts 
with considering the effects of the color octet, the QCD renormalization and the spin-orbit suppression. The dotted 
and dashed lines denote the results without considering the color octet and the QCD renormalization respectively. 
The experimental phase shifts [29, 30] are shown by black squares with error bars. 

Fig. 4: The theoretical predictions for the KN S, P and D-wave phase shifts. 

Fig. 5: The approximate expressions of the integral f(x) in Eq.(32) given by different values of the parameter 7 
which are plotted with the dotted lines. The real values of the integral is represented by the solid line. 

Fig. 6: Illustration of the effect of the spin-orbit suppression on the P-wave phase shifts. The solid lines represent 
the final results given by taking the parameter 7 = 0.45. The dotted lines denote the results given by 7 = 0.30. 

Fig. 7: The QCD effective coupling constants obtained from the one-loop renormalization. The solid, dashed and 
dotted lines represent the results given by the time-like momentum subtraction, the space-like momentum subtraction 
and the minimal subtraction respectively. 

Fig. 8: The QCD effective quark masses obtained from the one-loop renormalization. The solid line represents the 
result given by the time-like momentum subtraction. The dashed line shows the real part of the effective quark mass 
given in the space-like momentum subtraction. 



20 




P/ab(MeV) 



20 




200 400 600 800 1000 1200 

P/ab(MeV) 




50 100 150 200 

P /afo (MeV) 




P /afo (MeV) 



10 


i i i | i i i 

: (b) 


i i i i | i i i | i i i 


_ 5 

CO 
O 

B 

CO 




p ~ 

.... - r 13 : 






-5 




P T1 i 


-10 


■ . . . 


■ — i 

i . . . i , , , i . . N 



I I I I I I I I I I I I I I I I I I I I 3 

200 400 600 800 1000 



P /afo (MeV) 



